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ABSTRACT 

As part of the Panchromatic Hubble Andromeda Treasury (PHAT) multi-cycle program, we ob- 
served a 12' x 6.5' area of the bulge of M31 with the WFC3/UVIS filters F27bW and F336IV. From 
these data we have assembled a sample of ^4000 UV-bright, old stars, vastly larger than previ- 
ously available. We use updated Padova stellar evolutionary tracks to classify these hot stars into 
three classes: Post-AGB stars (P-AGB), Post-Early AGB (PE-AGB) stars and AGB-manque stars. 
P-AGB stars are the end result of the asymptotic giant branch (AGB) phase and are expected in a 
wide range of stellar populations, whereas PE-AGB and AGB-manque (together referred to as the 
hot post-horizontal branch; HP-HB) stars are the result of insufficient envelope masses to allow a full 
AGB phase, and are expected to be particularly prominent at high helium or a abundances when the 
mass loss on the RGB is high. Our data support previous claims that most UV-bright sources in the 
bulge are likely hot (extreme) horizontal branch stars (EHB) and their progeny. We construct the first 
radial profiles of these stellar populations, and show that they are highly centrally concentrated, even 
more so than the integrated UV or optical light. However, we find that this UV-bright population 
does not dominate the total UV luminosity at any radius, as we are detecting only the progeny of the 
EHB stars that are the likely source of the UVX. We calculate that only a few percent of MS stars in 
the central bulge can have gone through the HP-HB phase and that this percentage decreases strongly 
with distance from the center. We also find that the surface density of hot UV-bright stars has the 
same radial variation as that of low-mass X-ray binaries. We discuss age, metallicity, and abundance 
variations as possible explanations for the observed radial variation in the UV-bright population. 

Subject headings: galaxies: evolution-galaxies: individual (M31)-galaxies: stellar content-stars: 
evolution stars: horizontal branch 
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Many elliptical and large spiral bulges show enhanced 
ultra-violet (UV) flux towards their centers. This ex- 
cess from ~ 2000A to the Lyman limit is referred to as 
either the UV upturn, UV -rising branc h, or UV-excess 
(UVX). First discovered by 
the UVX converged on stel 



Code (1969), explanations of 
!ar phenom ena as higher res- 
olution telescopes came online j Bertola et al. 11995 ; King 



1995 



et al. || 1992 
spectra fro m HU 
et^p97 |. 



Brown et al. |1998| |2000| and fax^UV 
and FUSE became available (Brown 



1'he U VX is generally believed to be due to core helium 
burning stars (commonly known as horizontal branch 
stars; HB) and their descendants. These stars have 
evolved off the zero-age horizontal branch (ZAHB) and 
are following some combination of three evolutionary 
channels to becoming white dwarfs, depending on the 
fraction of their stellar envelope lost while they were red 
giants. Stars with a modest amount of mass-loss on the 
red giant branch (RGB) will leave the red side of the 
HB, lose a large convective envelope as canonical asymp- 
totic giant branch (AGB) stars, and become Post-AGB 
(P-AGB) stars. Intermediate amounts of mass loss will 
leave stars with smaller convective envelopes (which will 
be subsequently lost) and hotter effective temperatures 
than canonical AGBs. These stars will leave the AGB 
track early, becoming post-early AGB (PE-AGB) star^j 

18 This same effect can also arise in certain scenarios where stars 
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Stars that lose nearly all their envelopes on the RGB 
become extreme (blue) horizontal branch (EHB) stars. 
These stars will have envelopes which are too small and 
temperatures which are too high to reach the canoni- 
cal AGB line, and instead will beco me AGB-manque or 



post-EHB stars (for a full review see Greggio & Renzini 



1999| |Q'ConnelIfl999| ) 

Examples of these stars' evolutionary tracks and the 
hot side of the ZAHB arc illustrated in Figure [l]for both 
a Hertzsprung-Russell diagram and UV color-magnitude 
diagram (CMD). AGB-manque and PE-AGB stars col- 
lectively are long-lived, hot post-HB stars (hereafter re- 
ferred to as HP-HB). These HP-HB are 2-4 magni- 
tudes brighter in ^33614^ than EHB stars with similar 
F275W - F336W color. HP-HB stars have bolomet- 
ric luminosities around 10 2 ~ 3 Lq and have typical ef- 
fective temperatures of 15,000-20,000 K, compared to 
L ~ 3 L & and T c g = 10 4 K for HB stars (Bressan et al. 
2012 1. The initial masses of stars that become HP-HB is 



never more than 2 Mq, corresponding to main sequence 
( MS) lifetimes > 0.6 G yr. 
Brown et al| ( |1998[) (hereafter [B98| obtained HST 



Faint Object Camera (FOC) data in the M31 bulge, pro 
viding one of the most detailed discussions of these differ- 
ent types of UV-bright stars. However, these data only 
covered a small area (14" x 14"; see magenta box in Fig- 
ure^. Since radial gradients in age and metallicity are 
common in ellipticals and the bulges of large spiral galax- 
ies (e.g., |Carollo et al.||1993| |Davies et aL]|1993| [Trager 
et al |2000l|Peleiier ei a" liyOollGorgas et al'|iyy0lfl4iom-| 
sen fc Baum|1987 |Wirth|1981 1, we should expect radial 
gradients in the UV source population as well. These 
gradients can be used to constrain how age and metal- 
licity affect the stellar evolution leading to UV HP-HB 
stars. 

M31 is an excellent target for the study of radial trends 
in galaxy bulges. It hosts the most massive and metal 
rich bulge (e.g., Ferguson & Davidsen 1993) that can 
be resolved into incTivid ual U V so urces (as demonstrated 
by |Bertola et al. 1995 and B98), has a measured UVX 



( |Burstein et al.|1988 ), and estimated age-m etallicity gra- 



dient s derived from integrated line indices (|Saglia et al. 
20101 



In this paper, we present observations and analysis of 
the UV component of the bulge of M31, the first data im- 
aged in the UV by the Pan chromatic Hubble Andr omeda 
Treasury (PHAT) survey ( jDalcanton et al.||2012[ ). Sec- 
tion [2] presents our observations and their comparison to 
those of |B98| We investigate the radial properties of the 
UV sources and present a first broad comparison to new 
stellar models in Section [3] In Section |4j we discuss our 
findings in light of the possible causes of the UVX. Clos- 
ing remarks are presented in Section [5j All magnitudes 
quoted throughout this paper are in the VEGAMAG sys- 
tem. 

2. THE DATA 

2.1. PHAT Observations, Resolved Star Photometry, 
Astrometry 

As part of the PHAT program (GO-12058), we ob- 
tained UV imaging over a 12' x 6.5' region (2.6 x 1.4 kpc) 

burn their envelopes from the bottom up on the ZAHB, also leading 
to PE-AGBs. 




log T elI (A') 
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Fig. 1. — Sample stellar evolutionary tracks | |Bressan et al.|2012| l 
showing hot HB stars and their progenies' placement on tih, dia- 
gram (left) and CM D in UVIS filters (right) as seen at the distance 
of M31 (see Section |3.1| . Each evolutionary channel is illustrated 
and labeled. A canonical AGB track (black) is not shown on right 
pane l as it is fainter than F336W = 28. The P-AGB track is taken 
from [Vassiliadis fc Wood] ( |1994| l H-burning tracks for a mass of 
0.59Y Mq and Z.=U.Ult>. Dashed box shows £ 
other CMDs in this work. 



approximate range of 



around the center of M31 with HST/WFC3-UVIS. The 
UVIS channel has a pixel size 15/im and plate scale of 
0.04 arcsec/pixel. Images were taken in the F275W and 
F336W filters. Although the filters were chosen primar- 
ily for the study of massive stars in the star-forming disk, 
they also allow an excellent sampling of the UV-bright 
sources across the bulge. The area covered for this study 
is illustrated in Figure [2] and consists of a 3 x 6 grid of 
pointings with a 180° orientation flip between the two 
3x3 subgrids. The UVIS fields overlap by - 45". Two 
exposures were taken at each position in each filter, with 
a 1.9" dither between exposures to cover the UVIS chip 
gap. Total exposure times are 1010 seconds and 1350 sec- 
onds in F275W and i^336IU respectively. Fo r complete 
details of th e PHAT observing strategy see Dalcanton 
eFaL]([2012). 

We performed point spread function photometry on 
all UVIS pipeline processed data (.fit files) in the re- 
gion of interest using the software package DOLPHOTpj 
DOLPHOT is a modified version of HSTPhot ( polphm 
20001 that has been updated to include a specialized 
WFC3 module. Cosmic rays were rejected from the raw 
images using the IDL package lacosmic, which masks 
cosmic rays based on their very sharp edges. Further- 
more, our sharpness cuts (see below) removed most other 
cosmic rays from the photometry catalog. 

We have merged the independent photometric cat- 
alog for each of the 18 pointings into a single cata- 
log for the entire bulge. We have removed duplicate 
sources in overlapping regions as follows. First, given 
DOLPHOT's source positions and initial astrometric so- 
lution, we searched the overlapping regions for pairs of 
stars that are close in celestial coordinates. We used a 
search radius of 2 arcseconds which is large enough to 
compensate for HST pointing errors. We then produced 
a histogram of (ARA, ADec) vectors, which yielded a 
large but nearly uniform background of false matches 
plus a bump of correct matches at the RA, Dec offset 
between the two pointings. We measured these offsets 
for all pairs of pointings, then performed a least-squares 
fit to find the offsets and affine corrections (scales, rota- 
tions, and shears) between the pointings. By applying 

19 http:/ /purcell. as. arizona.edu/dolphot 
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Fig. 2. — Left: HST F336W (blue), F475W (green), and F8UW (red) composite mosaic of M31 bulge. Elliptical contours show 
locations of isophotes used as analysis region boundaries. The color scheme is used in other figures in this work and the median r adiu s 
of each isophote bin is superimposed. Data withi n the innermost (black dotted) region are excluded from the analysis (see Section |2.3| . 
The magenta box sh ows the footprint of the |B98| FOC imaging. Right: GALEX FUV (blue) and NUV (green) composite mosaic of M3l 
parmby et al.|[2006[ |. The red ellipse denotes the outermost (black) contour displayed in left panel. The footprint of the M31 disk field 
plotted in Figure [5] is shown as a red rectangle in right panel. 



TABLE 1 

Photometry Measurements 



Region 


Region 


Region 


Region 


median 


median 


-^stars 


90% 


90% 


Frac 


Frac 


area 


^min 


^max 


a 


a 




Comp. 


Comp. 


Cont. 


Cont. 




(arcsec 2 ) 


(arcsec) 


(arcsec) 


(arcsec) 


(pc) 




F275W 


F336W 


F275W 


F336W 





76 


0.0 


4.3 


2.1 


8.2 












1 


600 


4.3 


13.6 


8.9 


34.0 


280 


24.3 


23.9 


0.27 


0.22 


2 


3197 


13.6 


33.0 


23.3 


88.5 


607 


24.4 


24.0 


0.12 


0.14 


3 


9816 


33.0 


63.0 


48.0 


182.3 


856 


24.5 


24.4 


0.04 


0.09 


1 


34422 


63.0 


123.0 


93.0 


353.2 


1187 


24.5 


24.7 


0.02 


0.03 


5 


85143 


123.0 


214.0 


168.5 


640.0 


1404 


24.6 


25.1 


0.01 


0.01 



Note. — Region area is of each logarithmically spaced isophote; a is the semi-major axis (see Section |2.3ty , given 
in minimum, maximum, and median. The median radius is ca lculated directly from the isophotal boundary, assuming 
distance modulus (m — M)q — 24.47 l |McCoimachie et al.|2005| . -/V s t ars , are the numbers of UV-sources in each annulus 
(detected in both filters and brighter than the 90% completeness magnitude of the innermost bin); 90% Comp., percent 
completeness magnitude i n ea ch annulus. Frac Cont., fractional contamination at the 90% completeness magnitude of 
each annulus (sec Scction|3.4[. 



these corrections, we could put the sources detected in 
the separate pointings on a common local astrometric 
system. We then searched for sources that are measured 
in multiple pointings and merged their photometric mea- 
surements so that each star was represented by a single 
catalog entry. The local astrometric system was then 
tied to a global astrometric system by first aligning the 
WFC3/UVIS sources with our ACS/WFC sources, which 
are in turn aligned to a local reference catalog produced 
from CFHT imaging, which is in turn tied to 2MASS 
( |Skrutskie et~alp006| . 

The photometry output was then filtered to only allow 
objects classified as stars with DOLPHOT parameters 



signal-to-noise > 4, sharp < .075, crowd < 0.5, 
and round < 1.5 in both filters (see Dolphin||2000 for a 
detailed description of each parameter). 

To measure the completeness function, we performed 
100,000 artificial star tests (ASTs) for each observed 
field, and applied the above photometric cuts to the 
results. We find the data to be at least 90% com- 
plete in each radial bin down to F275W = 24.3 and 
F336W = 23.9. We limit our analysis to this innermost 
90% completeness limit to ensure that completeness is 
at least 90% at all radii. Completeness magnitudes are 
listed in Table Q] for each radial bin. We made one fur- 
ther cut to limit the contamination of hot Milky Way 
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Fig. 3. — A comparison of PHAT UV and FOC spatial over- 
lap. Left: UV CMD for the [598] dataset. The FOC sources 
with PHAT counterparts are shown in red (filled circles), and FOC 
sources without PHAT counterparts are open circles. The PHAT 
data appear incomplete for the bluest sources with FnbWpoc ~ 
F275W FO c < -0.75. Right: UV CMD for t he P HAT dataset 
extracted from the smaller FOC footprint of B98 Stars that 
have matched counterparts in the |B98| catalog are plotted in red 
(filled circles), and unmatched PHAT sources are open circles. The 
completeness limit of the P HAT data is roughly 0.3 magnitudes 
brighter in F275W than the B98 FOC observations. The dashed 
lines correspond to|B98 detection limit. 

(MW) foreground stars by removing bright red stars with 
F336W < 19 and ^2751^ - .F336W > 0.8 (see Section 
2.31. 



2.2. Comparison with FOC data 

As a c onsist ency check, we compared our photometry 
to that of |B98| which carried out a comparable UV stellar 
population anal ysis, b ut over a much smaller area than 
presented here. |B98 s imaging was taken by the FOC 
in its F175W and F275W filters, which are suited for 
the detection of stars with effective temperatures up to 
T ff ~ 40, 000 K. To compare these datasets, we began by 
correcting the FOC photometry to be 0.5 mag br i ghter 



Brown et al.| ( |2000 1 
We made no conver- 



as recommended by Section 4.4.1 in 
to account for calibration revisions 
sion between the FOC and WFC3 filter systems, beyond 
transforming the |B98| data from its tabulated STMAG 
values into the VEGAMAG system. We then trans- 
formed the |B98| catalog positions onto our astrometric 
system using linear shifts determined by matching bright 
sources (F275W < 22) in the similar F275W passbands. 
Once aligned, we constructed comparison catalogs with 
identical spatial coverage defined by the 14" x 14" FOC 
field of view (see Figure [2]) and excluded a circular 2.5" 
radius in the center to exclude the most crowded regions. 
Finally, we performed cross-catalog matching using a ~ 2 
pixel (~ 0.07 arcsec) search radius, and accepted matches 
with F275W magnitude differences < 0.7 mag. If mul- 
tiple sources fulfilled the match criteria, we chose the 
match with the smallest magnitude difference. 

The results are presented in Figure [3] We were able to 
match 65% of the B98 catalog to WFC3 sources m our 
catalog. A majority of the unmatched |B98| sources fall 
in the fainter (F275W FO c > 24.5), bhierjF175W FO c ~ 
F275Wfoc < -0.5) portion of the FOC CMD, showing 
the |B98| dataset 's slightly shorter wavelength sensitivity. 
Overall, however, this comparison showed that these two 
datasets probe comparable stellar populations, with the 
PHAT dataset providing superior spatial coverage, and 
the FOC sample providing a better probe of the very 
hottest stars. 



2.3. Radial Binning of UV Sources & Possible 
Contaminants 

We divided the data into bins to analyze the properties 
of the UV sources as a function of radius. These bins 
were defined using logarithmically space d isophotes, as 



meas ured on a 3.6/im image from Spitzer (Barmby et al 
20061. The innermost contour (semi-major axis distance 
of ~ 4") was set to exclude the nuclear region, due to 
the relatively lower quality of our photometry, while the 
outermost contour (214") was set to avoid major dust 
lanes and regions of recent star formation that lie along 
the inner arms/ring in the disk. Table Jl] lists the proper- 
ties of the bins, which are illustrated in Figure [2] Their 
CMDs are shown in Figure [4] 

Because metal poor MS stars can occupy the same re- 
gion of CMD space as UV-bright stars, we needed to as- 
sess potential effects of contamination on our study. We 
examined the contamination from the star-forming disk 
of M31 by comparing the radial gradients of stars in two 
regions of the CMD. The first CMD region was chosen to 
occupy the same color-magnitude space as a MS domi- 
nated disk field (galactocentric radius ~ 6.5 kpc), shown 
as a red rectangle in the right hand image in Figure 
(22 < F336W < 23, -0.13 < F275W - F336W < 0.28 
The second CMD region was selected to be too blue to 
have contamination from MS stars (22 < F336W < 
23,-0.35 < ^275^ - F336W < -0.50). The stellar 
gradients of these regions showed a nearly constant ratio 
within uncertainties, indicating MS contamination from 
the M31 bulge is small across the analysis region. Fur- 
thermore, bulge-disk decomposition suggest that while 
the fraction of disk light increases five-fold between the 
inner and outer annuli, the disk never contributes more 
than 12.5%, even in the outer most analysis region (How- 
ley et al, in prep). 

There is also a possibility that some of the redder stars, 
(F275W - F336W > 0.8; see Figure [§ are hot MW 
foreground stars. To esti mate the MW conta mination we 
have made TRILEGAL ( Girardi et al.|2005 ) simulations 
of the MW for the field ot view contained, in the entire 
analysis region. Most of the simulated foreground stars 
are redder than F275W - F336W = 0.8 and brighter 
than F336W = 19. We excluded 69 sources brighter 
and redder than these limits from our analysis. Bluer 
and fainter than these limits, our simulations suggest a 
negligible contribution (~ 40 stars) of foreground stars 
with F33QW > 19. 



2.4. Integrated Light and Ancillary Imaging 

As a way of assessing bulk properties of the defined 
analysis regions, we calculated total integrated fluxes for 
each analysis region from the F275W and F336W imag- 
ing, as well as from the Spitzer/IRAC 3.6^tm image. For 
the HST imaging, we first assembled an aligned mosaic 
of the UV imag ing by means of the M ultidrizzle task 



within PyRAF (Koekemoer et al. 2002 1 using the as 



trometry solution obtained from the stellar catalogs and 
with sky-subtraction disabled. Once we obtained this 
combined mosaic, we performed aperture photometry of 
the regions on the combined image. To estimate the sky 
level for the F275W and F336W mosaics, we performed 
photometry on outer disk images (galactocentric radius 
of ~ 15 kpc) obtained with identical observing parame- 
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Fig. 4. — UV Color-magnitude diagrams of each analysis region. Dashed lines mark the 90% completeness limits and solid lines mark 
the 50% completeness limits (see Section ^J. We limit out analysis to sources brighter than the 90% completeness limits in the innermost 
(green) region. The median projected ramus of each distance-corrected logarithmically-spaced isophote is labeled in the upper left. The 
color scheme follows Figured 
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Fig. 5. — Nature of the bright UV point sources in the central 
bulge (left) vs. star-forming regions in the field (right). The two 
panels on the left show the UV CMD for stars within 120 pc of the 
center; the two panels on the right show CMDs of the disk field 
(galactocentric radius ~ 6.5 kpc; see red rectangle in the right panel 
of Figure [M . These CMDs are compared to MS tracks for young 
massive stars in the top panels; they are compared to P-AGB and 
HP-HB tracks for evolved stars in the bottom panels. It is appar- 
ent that the (reddened) MS tracks ar e a g ood match for the disk 
CMDs, but not the bulge (see Section |3.1[ l. Markers are placed on 
the three most massive HP-HB tracks at 10 4 yr intervals. Two red- 
dening vectors with Ry =3.1 and 5 are also shown for 1 magnitude 
extin ction in F336W. MS tracks assume Ay = 0.99 l |Kang et al.| 
|2009| ) and bulg e tracks assume a pu re MW foreground extinction 
oF7TV = 0.206 | |Schlegel et al.|1998| l. 

ters as part of the PHAT survey and confirmed that the 
true sky background level is ~ at these very blue wave- 
lengths. As such, there was no need for sky subtraction 
as part of our aperture photometry. For the Spitzer im 



age, we used available mosaics (see Barmby et al. 2006 
for details) and similarly performed aperture photometry 
for the analysis regions. 
Table [2] includes the su rface density of HP-HB and 



P-AGB stars (see Section 3.3 for a discussion on how 



the populations are separated) , their resolved UV flux, 
the total resolved UV flux, and total integrated UV 



flux in each filter, in each annulus. Also tabulated is 
the total integrated flux within each annulus obtained 
from Spitz er/IRAC 3.6/im imaging (F 36 ^ m , Barmby 
et al.|2006| and the bulge-to-disk ratio {B/D) integrated 
within the radial bin using a 2D, two-component model 
of an exponential disk and Sersic bulge (Howley et al. , in 
prep). 

3. ANALYSIS 

3.1. Radial Color-Magnitude Diagrams 

Figure [5] shows a UV CMD of the inner two radial bins 
(15-120 pc; left panels) and a disk field (right panels, 
red rectangle in Figure [2]) composed of the entire Brick 9 
(galactocentric radius ~ 6.5 kpc) from the larger PHAT 
dataset. We superimposed stellar evolutionary tracks 
converted to the WFC3/UVIS p hotometric system fol- 
lowing the Girardi et al. (2008) bolometric corrections 
updated with the latest VVF(J3/UVIS filter throughputs. 
We assumed a distance mo dulus of (m — M) = 24.47 



( |McConnachie et al.||2005[ ). When plotting young MS 
tracks, we adopted an interstellar extinction of Ay — 
0.99 f or the disk field, ba sed on the median E(B-V) value 



(2009) for star forming regions in the 
. Although we are not currently able 



from |Kang et al. 
disk analysis held 
to derive the foreground interste llar extinction for th e 
bulge, we adopted the value from Schlegel et al. (1998), 
Ay — 0.206. These Ay values were converted into the 
F275W and F336W extinction values using the coeffi- 
cients Ap275 w/Ay = 2.05 and A fmg w/Av = 1.67 de- 



rived from the |Cardelli et al.| ( |1989[ ) and |0'Donnell| ( |1994j 

(shortward of U -band) extinction curve with Ry = 3.1 
applied to a T e g = 15, 000 K star. 
The evolutiona r y trac ks come from two sources: 



Vas- 



siliadis fc Wood| (n994) for H-burning P-AGB starTof 
masses > 0.6 MQWithZ = 0.016 and Y = 0.25 and the 
most recent Padova library ( Bressan et al. |2012 ) for the 
HP-HB starj^°] For the latter, the chemical composition 
is cither Z = 0.07 and Y = 0.389 with an a-enhanced 



composition typical of bulges (adapted from Bensby et al. 
20T0| ) for the bulge field (left panels), or Z = 0TD2" and 
y = 0.285 with solar-scaled composition for the disk field 
(right panels). The metallicity chosen for the bulge is 
an extreme case to allow for high helium content. It 

20 P-AGB stars with masses M > O.9A/0 are not included in 
Figure [5] because their progenitors are too massive and short-lived, 
and thus not expected in the bulge. 
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may be possible to populate the EHB and its progeny 
in similar CMD space with evolutionary tracks of lower 
metallicity by increasing mass loss rates. The detailed 
analysis required to examine this possibility is forthcom- 
ing in Rosenfield et al. (in prep). The calculations for 
the tracks plotted in Figure [5] include all kinds of single 
stars that can give rise to UV-bright stars, namely the 
younger MS, the intermediate-age Post-AGBs, and the 
old EHB and HP-HB stars. 

It is evident from the left panels of Figure [5] that the 
CMD of the inner 15-120 pc of the bulge shows quite dif- 
ferent characteristics than the disk CMD in the opposing 
panels. The disk field shows a well populated diagonal 
sequence that is mostly redder than F275W — F336W = 
—0.3 and that becomes systematically redder at fainter 
magnitudes. In contrast, the bulge contains a significant 
population of stars bluer than F275W - F336W = -0.3 
that does not appear to be an under-populated version 
of the disk field. 

Comparing to models, the stars in the young, star- 
forming disk match the solar MS evolutionary tracks, as 
expected. However, the higher metallicity MS tracks fail 
to match the bulge population - with its wide spread 
in the CMD for F336W < 23, and its remarkable con- 
centration of stars at F336W > 23 - indicating that 
the UV-bright stars in the bulge are not primarily MS 
stars. In contrast, the HP-HB stellar evolutionary tracks 
provide a good match to the bulge, but fail to match the 
disk population. This match suggests that the UV-bright 
stars in the bulge are primarily a mix of PE-AGB and 
fainter AGB-manque stars (red tracks in Figure [5] red 
and blue tracks, respectively, in Figure [I]) . While more 
luminous P-AGB stars (purple tracks in Figures [I] and 
[5]) with masses > 0.6Mq could potentially be detected, 
m practice their evolution is sufficiently fast that they 
make no significant contribution to the UV CMD. 

If we examine the CMDs across the entire bulge (Fig- 
ure El), there is a clear radial trend. The number of stars 
brignter than F336W ~ 23 decreases towards the center 
as does the area of each isophotal region. Dimmer stars, 
however, do not show this trend, and in fact, their num- 
bers remain high at all radii. Both of these trends are 
probably driven by changes in the relative proportions of 
P-AGB and HP-HB stars, which in turn depends on the 
changes in the progenitor stars as a function of radius 
(though the P-AGB do not contribute to the faint gra- 
dient). We analyze these trends more extensively below. 

3.2. Radial Luminosity Functions 
Figure [6] shows the luminosity functions (LF) for 



M/L ratio of 1.09 derived using t he color-dependent rela- 
tions from Bell & de Jong ( 2001[), and 3 .6/nm to if -band 



contaminant-corrected (see Section 2.3) resolved UV 



sources within our adopted annuli. We removed the ef- 
fects of changing radial stellar density by scaling the LF 
by each region's stellar mass. We estim ated the stellar 
masse s from the integrated 3.6/mi flux (|Barmby et al. 



2006 1 , which should scale linearly with stellar mass, up to 
a multiplicative cons tant p^] Following the procedure of 



Barmby et al. fl2006|), we assumed an optical bulge color 
oiB-K = 1.8 ( [Walterbos & Kennicutt|1987| , a iVband 



21 Isochrones from [Gfrardi et al.| | |20io| and ref. therein) con- 
firm that the expected 3.b/im mass-to-hght ratio does not change 
significantly for stellar populations of ages older than ~ 4 Gyr and 
metallicities Z> 0.001. 



M/L ratio conversion from Oh et 
in an 3.6/im M/L ratio of 0.95. 



al. 



In 



(2008K resulting 
jolwe made 
sak 



? igure 



completeness corrections to see if each LF peak was not 
simply a manifestation of incompleteness of the data. To 
do so, we removed our magnitude limit criteria (see Sec- 
tion [2| and divided the sources of good quality detected 
in both filters by the completeness fraction found from 
the ASTs in 0.2 magnitude bins. 

The LFs look similar at bright magnitudes for all radii, 
but show significant radial variations at fainter magni- 
tudes. First, we see a dramatic trend at faint magni- 
tudes where faint HP-HB stars (F275VF, F336IV > 23) 
become increasingly abundant in the inner bulge, relative 
to the overall stellar mass. Second, we see a weaker trend 
for stars with magnitudes between 21 — 23, which also 
show increasing contributions toward the inner bulge. 
These are likely P-AGB stars of a continuous distribu- 
tion of masses less than ~ 0.6Af Q , suggesting that the 
outer regions are under-abundant in P-AGB stars. 

To help interpret these results, we overlay a syn- 
thetic LF (dotted in Figure [6]) made of low-mass Z = 
0.07, Y = 0.389 tracks from t he Padova stellar evolution 
library (Bressan et al. 2012) evolved from the ZAHB, 
assuming that the mass distribution of stars along the 
ZAHB is flat for all ZAHB masses below 1.0 M e . This 
approximation is reasonable since the range of ZAHB 
masses that actually contributes to the LF is very nar- 
row (~ 0.45 - O.58M ). 

The simulated helium-burning LF reproduces the 
shape and color of the bump in the LF fainter than 
22.5 mag, especially for the inner M31 radial bin. This 
bump can therefore be interpreted as being made by the 
HP-HBs. Deeper observations would likely reveal a much 
larger bump in the LF due to the EHB. 

In addition, the flux density from the detected stars 
does not vary with magnitude except for the inner re- 
gion. Essentially all stars between 20 and 24 in both 
F275W, F336W contribute equally to the flux density 
with the exception perhaps of region 1, where the con- 
tribution tends to increase towards fainter magnitudes. 

3.3. Luminosity Contribution of HP-HB and P-AGB 

Although P-AGB and HP-HB stars are individually 
bright in our UV bandpasses, together they only make 
up a small fraction of the integrated UV light (< 2% 
in F275W). A vast majority of the UV light must be 
emitted by fainter populations, such as EHB stars that 
remain undetected in our resolved photometry. Theoret- 
ical calculations (Rosenfield et al., in prep.) that match 
the resolved star populations predict a flux contribution 
from faint unresolved stars that is comparable to mea- 
sured flux in unresolved stars. 

Since the P-AGB and HP-HB stars we detect are a very 
small fraction of the UV light (see Table [2]), the total in- 
tegrated UV flux must be coming from some other pop- 
ulation of fainter, but more numerous stars. Our model 
LF suggests a large population of EHB stars fainter than 
mag ~ 25. These stars are likely the main source of 
the UVX (c.f., B98). These stars could be detected with 
deeper UVIS observations and with bluer filters, provided 
that the crowding limit was sufficiently faint to detect 
these numerous stars. 
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Fig. 6.— Number of resolved HP-HB stars (left panel: F275W; 
right panel: F336W) as a function of magnitude, scaled by each 
regi on's stellar mass found by the 3.6/xm M/L ratio (see Section 
|3.2| l, in bins of 0.2 magnitudes. The observed LFs are corrected for 
completeness and shown with thin lines fainter than the 90% com- 
pleteness magnitude in the respective annulus (i.e., a correction by 
more than a factor of 1.1). The color scheme is the same as Fig- 
ure [2] with green indicating the innermost bin. Long lived HP-HB 
stars are shown in model LF (dotted) evolved from the ZAHB with 
a flat distribution of mass below l.OM , and the Z=0.07 Y=0.389 
tracks (cf. Figure [5}. The model is normalized to the innermost 
region's peak LF m F275W. The over-abundance of stars with 
F275W > 23 in the innermost regions of M31 are associated with 
HP-HBs. The bright LF tails are due to a seemingly flat distribu- 
tion of P-AGB stars (see Section |3.5| . The increase in the model 
LF at faint magnitudes (> 24.5) is due to numerous HB stars, of 
which the hottest likely dominate the UVX. 

To further quantify the radial gradients of stellar pop- 
ulations, we divided UV-bright stars into HP-HB and 
P-AGB stars. We separated the HP-HB stars from 
P-AGB stars using the highest mass evolutionary track 
that produces a HP-HB star at Z = 0.07 and Y = 0.389 
(M = O.58M ; see the brightest red track in the bot- 
tom left panel of Figure [5]) for the Padova stellar evolu- 
tion library. More massive evolutionary tracks than this 
critical mass produce stars with enough envelope mass 
to become canonical AGB stars; this behavior has little 
metallicity dependence, and does not app ear to change 
with He-abundance from 0.26 < Y < 0.46 (Bressan et al. 
2012|. We have classified stars brighter than this track 
as P-AGB. As a faint limit of UV-bright stars, we use 
our innermost region 90% completeness magnitude (see 
Section pi}. 

The left panel of Figure [71 shows the surface brightness 
profiles of the HP-HB, P-AGB, and integrated UV light. 
In both filters, the surface brightness increases toward 
the center for both HP-HB and P-AGB. The right panel 
of Figure [JJ shows the ratio of resolved UV stars, and 
indicates that the flux from HP-HB stars increases with 
respect to the P-AGB 



A best fit Sersic profile from Courteau _et al. (20111 
is overplotted on the left panel of Figure 7[ with arbi- 
trary n ormalization in surface brightness. Courteau et al. 
pOll) fit UVBRI, 2MASS, and Spitzer IRAC surface 
brightness profiles of M31 to find the Sersic bulge index 
of n = 2.2 ± 0.3 and effective radius of R e = 1.0 ± 0.2 
kpc. We did a simplistic least-squares fit to a Sersic pro- 
file using the surface brightness of our total integrated 
light in F336W imagery to find He F33ew = 21.7 mag/// 2 , 
n F336W = 2.4, and R eF336W = 2.1 kpc. These numbers 
are consistent with the steeper scaled gradient in UV 
sources. 

3.4. Possible Contributions from Stellar Blends 

The bulge of M31 has a strong surface brightness gra- 
dient in the UV, which must be due to a dramatic in- 
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Fig. 7. — The radial distribution of UV sources. Left Panel: 
Surface brightness as a function of radius of F275W (solid, black), 
F336W (solid, red), HP-H B (dotted), and P-AGB (dashed). The 
best fit Sersic profile from |Oourteau et ah] ( |2011^ is overplotted 
with £t e shifted to F336W='2'6 (gray). Kight Panel: Resolved flux 
ratio of HP-HB to P-AGB stars in each filter for each annulus. 
Shaded regions show Poisson uncertainty in number counts. In 
both panels, red denotes the -F336W filter and black denotes the 
F275W filt er. The characterist ic metallicity as a function of radius 
reported by Saglia et al. ( 20ll)]l is indicated on the top axes. Both 
the HP-HB and J-'-AlitS populations increase toward the nucleus of 
M31, but the HP-HB stars show a dramatic increase in flux with 
respect to the P-AGB stars in the innermost regions. 

crease in the number density of stars towards the center 
of the galaxy. As we show in Table [2j the majority of 
these stars are unresolved in our observations, but they 
could possibly lead to the detection of spurious sources, 
if sufficient numbers of the fainter stars were blended 
together to rise above our detection threshold. These 
spurious sources would dominate at fainter magnitudes 
and would be expected to show a strong radial gradient 
due to the increased crowding in the center, and thus 
could potentially mimic our observed gradient. While 
we have minimized the likelihood of this contamination 
by excluding the most crowded regions of the bulge and 
by choosing a relatively bright limiting magnitude for 
our analysis (F336W = 23.9, .F275VF = 24.3), it is still 
worthwhile to confirm that this issue is not affecting our 
conclusions. 

We have ruled out this possibility by considering two 
different cases. For the first, we assume that the luminos- 
ity function observed at bright magnitudes continues as 
a power-law to faint magnitudes. This assumption puts 
large numbers of stars just below the detection thresh- 
old, where the chances are maximized that a star will 
blend with an undetected source and rise above the de- 
tection threshold. We then simulate the observed lumi- 
nosity function as follows. In a series of magnitude bins, 
we randomly select stars from the same catalog of ar- 
tificial stars used to calculate our completeness limits. 
The number of stars drawn in each magnitude bin is de- 
termined by our assumed power-law luminosity function. 
We then use the recovered properties of the artificial stars 
(i.e., the magnitudes the stars were recovered with, and 
if they were detected at all) to generate the luminosity 
function that would be observed. Figure [8] shows the in- 
put and recovered luminosity function in each of the ra- 
dial bins, plotted as black and colored lines, respectively. 
As expected, blending leads there to be an artificial up- 
turn in the luminosity function as one approaches the de- 
tection limit, with the effect being more dramatic in the 
most crowded inner annulus. However, these effects only 
become appreciable below our adopted magnitude cuts 
(vertical dashed lines). At our adopted analysis limit of 
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Fig. 8. — Input (black) and recovered (blue: F275W; green: 
F336W) luminosity functions in each of the radial bins. Close 
to the detection limit, blending causes an artificial upturn in the 
luminosity function. The effect only becomes appreciable below 
our adopted magnitude cuts (dashed lines) and is never above 27% 
(see Table ^ for exact values in each radial bin). 

F336W = 23.9 and F275W = 24.3, the observed lumi- 
nosity function never has more than 27% contamination 
from blends (see TablefTlfor exact contamination levels in 
each bin) . Given that the luminosity function is observed 
to vary a factor of 100, it is unlikely that the observed 
radial trend is produced by radially-dependent blending 
from a power-law luminosity function. 

We have also considered a second, somewhat less likely 
source of blends. It is believed that the majority of the 
unresolved flux in our images comes not from sources just 
below our detection threshold, but from a large popula- 
tion of long-lived EHB stars (e.g., B98) that are more 
than a magnitude fainter than our detection threshold. 
If so, these sources must have a very high surface density, 
which could sometimes lead to multiple EHB stars blend- 
ing together within a single resolution element. A ssum- 
ing a typical EH B magnitude of F33QW = 26.5 (Bres- 
san et al. 2012), and a surface brightness of ~18 mag- 
nitudes arcsec -2 , we can calculate the mean number of 
EHB stars per square arcsecond. Adopting this mean, 
we can then calculate the Poisson probability that there 
will be a sufficiently large upward fluctuation in the local 
density (within the HST resolution element) to produce 
a spurious source at F336W = 24. We then multiply by 
the number of independent resolution elements to calcu- 
late the contamination in the inner most annulus. When 
we do so, we find that we expect no more than 0.1 spu- 
rious stars brightward of F336W — 24, in contrast to 
the 164 stars found between 23.5 and 24 magnitudes in 
region 1. We therefore rule out blends from EHB stars 
as a significant source of contaminants to the observed 
luminosity function. 

These two tests indicate that it is highly unlikely 
that false blended sources are significant contaminants 
brighter than the magnitude limit we adopted for our 
analysis. 

3.5. Comparison to Other Evolving Populations 

Figure [7] establishes that the resolved flux of HP-HB 
increases with respect to the P-AGB sources toward the 
nucleus. We now look in depth at the radial variations of 
these populations and expand the discussion to include 
other stellar populations. 

The HP-HB sources discussed here are not the first 
population of astrophysical sources to show enhance- 
ment towards the inner regions of M31 's bulge. In their 
study o f Chandra X-ray point sources, |Voss fc Gilfanov| 
( 2007b ) show an increase in the low-mass X-ray binary 
(TMXB) population at small radii within the bulge of 
M31. These authors conclude that this enhancement 
can be accounted for by a population of dynamically 
produced LMXBs, which explains why the distribution's 
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Fig. 9. — Surface density of U V-bright sources ( blue = HP-HB, 
red = P-AGB), LMXBs (green, fvoss fc Gi lfanov 2007b), binned 
according to analysis region with Poisson uncertainties. The 3.6/im 
luminosity (black dashed) is plotted in different units and follows 
the right axi s. The characteristi c metallicity as a function of radius 
reported by |Saglia et al,| ( |2010} is indicated on the top axis. The 
least-squares nt power law index of each stellar component is noted 
above its profile with its uncertainty in the fitting. HP-HB (blue) 
and P-AGB (red) stars follow a different surface density profile 
than that of the 3.6/^m light. However, the HP-HB stars seem to 
follow a profile consistent with the LMXBs. 

profile follows the square of the stellar density (pi). 

For comparison, we bi nned the X-ray sources from the 
Voss & Gilfanov ( 2007a I catalog according to our analy- 
sis regions, excluding those sources that are either iden- 
tified as non-LMXBs or that are associated with globular 
clusters (as their origin is different than that of unclus- 
tered LMXBs). The resulting number density distribu- 
tions for the LMXBs are shown in Figure [9] We find the 
intriguing result that the excess in UV-bright sources is 
consistent with the power law index of the LMXBs. 

We have also evaluated the radial trends in the density 
of PNe, as these stars are likely descendants from the 
P-AGB populations and might b e expected to s how a 



similar radial trend. We use the Merrett et al.| ([2006 1 
catalog of 2615 PNe in M31 detected by 5007A emission 
and find that the distribution of PNe across the bulge 
is roughly constant with radius, and does not follow the 
trends of HP-HB, P-AGB, or LMXBs. However, the 
Merrett catalog is not expected to be complete in the 
innermost regions. If the true radial distribution of PNe 
does follow the gradient in the P-AGB population, then 
the inner regions must be incomplete by a factor of ~ 5.5. 
We discuss the implications of these results in Section [4] 
below. 

3.6. Stellar Population Lifetimes Inferred from 
Observations 

Now that we have detected gradients in HP-HB stars, 
it is useful to make a rough estimate of what fraction 
of stars are expected to become U V-bright through this 
channel as a function of radius. IRenzini &: Buzzonil 



Bright UV Stars in the Bulge of M31 



TABLE 2 

Radial Gradients 



Region 


B/D 




Integrated Flux 


Resolved Flux 


HP-HB Resolved Flux 


Surface Density 






F275W 


F336W 


F275W 


F336W 


F275W 


F336W 


HP-HB 


P-AGB 






(Jy) 














(arcsec -2 ) 


(arcsec -2 ) 





40.7 


0.27 


















1 


40.7 


1.16 


25.26 


98.33 


0.96 


0.75 


0.46 


0.39 


0.383 


0.082 


2 


33.5 


3.74 


89.77 


351.08 


2.23 


1.85 


1.01 


0.84 


0.150 


0.040 


3 


20.7 


6.77 


181.40 


681.08 


4.01 


3.45 


1.51 


1.23 


0.062 


0.025 


1 


12.0 


12.99 


373.58 


1315.45 


5.70 


5.12 


2.14 


1.74 


0.024 


0.010 


5 


7.0 


16.71 


608.05 


1802.78 


6.96 


6.73 


2.52 


2.03 


0.011 


0.005 



Note. — First two columns are the same as in Table [T] followed by B/D, the bulge- to- disk ratio of luminosity within each 
annulus (Howley et a l., in prep); i^Ve/im, region integrate a flux within each annulus obtained from Spitzcr/IRAC 3.6/j-m imaging 
jBarmby et a~ Next arc integrated and resolved flux measurements in each UVIS filter, followed by the resolved flux of 



HP-HB stars. All UVIS fluxes are in units of 10" 
populations are tabulated (cf. Figure fol. 



( 1986 ) provided a powerful tool to relate the number of 
evolved stars in a population to their lifetimes, with only 
a few conservative assumptions. Using the fact that most 
fuel is consumed by stars fusing hydrogen or helium, and 
that the initial mass function is not too steep, they relate 
the number of stars, Nj; in any evolved state, j, to the 
average amount of time (in years) the star spends in that 
state (tj): 



N ]= L T B{t)(t 3 



(1) 



where Lt is the total bolometric luminosity of the stellar 
population in Lq and B(t) is the rate of stars leaving the 
MS in units of Lq^ 1 yr _1 . 

Equation [l] is strictly valid only for a stellar popula- 
tion of single age and metallicity that happens to contain 
HP-HB stars. Any subregion of the M31 bulge, instead, 
is likely to contain a range of populations, with some 
spread in their ages and metallicities. Since just a frac- 
tion / of these populations will likely produce HP-HB 
stars (e.g. the oldest, most metal rich, and/or of high- 
est helium content), equation [l] can be conveniently re- 
written as 

/ = ^hp-hb ^ 
LtB(£) (£hp— hb) 

where f = fT HF ~ HB,obB . This equation can be solved for 

■> JV H P-HB,ppcd ^ 

/, allowing us to solve for the fraction of evolving stars 
that pass through the HP-HB channel. 

We have previously measured A^np-HB.obsj as listed 
in Table [2| We also adopt a value of B(t) = 2.2 x 
lCP n LQ~ yr^ 1 , derived from solar metallicity evolution- 
ary models. This number sho uld be accurate to 10% fo r 
populations 10 Gyr and older ( |Renzini k. Buzzoni| 1986). 
For the average time in the observed state, £hp-hb, we 
calculated the mean CMD crossing time of the relevant 
Padova stellar evolution tracks shown in figure [5] to find 
(*hp-hb) - 2 x 10 6 yr for the HP-HB. 

To calculate Lt, we convert the 3.6/im flux (corrected 
to exclude any disk contamination using B/D; see Table 
v3 to bolometric luminosity. As mentioned above, 3.6/im 
mix traces the stellar population of the bulge, modulo a 
bolometric correction: 



Li 



L3.6fj.111 x 



^bol 



L 



3.6/iin 



(3) 



We calculated the bolometric correction by simulating 
the integrated light of simple stellar popultations over a 



Finally, the surface densities of the resolved UV-bright 



range of ages and metallicities using isochrones from Gi 



rardi et al.| ( 2010 an d ref erences therein). As mentioned 
m pas sing in Section 3^ isochrones from |Girardi et al. 
(20101 confirm that the expected 3.6/im flux of star light 
does not change significantly for stellar populations of 
ages older than ~ 4 Gyr and metallicities Z> 0.001. 

With the bolometric correction, BC = 2.60 and adopt - 
ing M 0jbol = 4.77 and M , 3 . 6Mm = 3.24 ( |Oh et al.|2008| , 
we converted to Lt, 



L 



bol 



L, 



10 



-O.4(M 0ibol -M 0i3 . 6 ^ m -BC) 



2.7 



3.6/iin 



(4) 



Combining equations [2j [3j and substituting the values 
above for B(t) and (£hp-hb), we solve for /, 



/ = 1.7 x 10' 



.JVi 



HP-HB, obs 



(5) 



Substituting in number s fro m Table [2l we plot / as a 
function of radius (Figure 10 1. We find that the fraction 
of stars that go through the HP-HB channel is small, 
but that the fraction inc rea ses systematically towards the 
center of M31. Figure [l0| shows that nearly 3% of the 
MS turnoff (MSTO) stars in the innermost regions of 
M31 become HP-HB stars. This implies that 97% of the 
MSTO stars become canonical P-AGB stars. In contrast, 
in the outermost regions of the bulge, fewer than 1% of 
MSTO stars become HP-HB stars. Note that these -3% 
of HP-HB stars appearing as bright resolved UV sources 
belong to the same population of the fainter EHB stars 
that, according to evolutionary models, do mak e a l arge 
fraction of the unresolved UV light (see Section 3.3). 

The fraction of stars that go though the HP-HB chan- 
nel in the inner most regions is consi stent with what is 

Our findings are 
who demon- 



found in M32 by Br own et al. (2008) 



also consistent with Brown et al. ~| 



997) 



strated that galaxies with moderate U VX have about 
2% of the main sequence population passing through the 
EHB. 

This test has shown that, even in the center of M31, 
very few stars from the bulge population are necessary 
to become HP-HB and explain the UV-bright stars we 
detect. What causes an RGB or EHB star to eventually 
become a HP-HB star only affects a small percentage of 
the underlying stellar population, but that small likeli- 
hood varies strongly with radius, for reasons we discuss 
below. 
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Fig. 10. — Population variations in the central ~ 0.5 kpc of M31. 
Shown (solid, blue) is the inferred fraction of evolving MS stars that 
must go through the HP-HB channel to explain the HP-HB number 
count profile, given the integrated light properties. This fraction 
rises by about a factor of 3 towards the center. Shaded regions 
denote Poisson uncertai nty. The Lick -Index based metallicity and 
its dispersion derived by |Saglia et al.| j2010| is shown on the right 
axis (dotted, black). 

4. DISCUSSION 



Our LFs and CMDs allow us to confirm IB98l s claim 
that HP-HB stars are the likely evolutionary channel for 
the production of most of the UV-bright stars in the M31 
bulge. However, ~ 20% of brighter UV stars in the inner 
two radial bins (15-120 pc) are compatible with P-AGB 
tracks (see Figure [9]), expected for stars with negligible 
mass loss on the KGB. 

We detect clear gradients in the number and types of 
UV-sources. We see a steep gradient at faint magnitudes 
(F275W, > 23) that likely arises from an in- 

crease in the HP-HB population towards the center of 
the bulge. These faint stars contribute about 60% of the 
resolved flux of the inner two radial bins in both F275W 
and i^336VU (see Table[2) and roughly 2% of evolving MS 
stars pass though this long-lived HP-HB channel before 
becoming white dwarfs (see Figure 10), the remainder 
pass though the short-lived P-AGB ' However, the de- 
tected stars only make up ^2% of the total flux in the 
F275W filter indicating that we have not yet detected 
the stars directly responsible for the UVX. These miss- 
ing stars are likely to be faint EHB stars, as suggested 
by|B98l 

We now outline possible causes of the detected radial 
gradients of HP-HB stars. 

4.1. Gradients in Stellar Ages 

The UVX is expected to increase with time for an aging 
stellar population. As a population ages, the mass of the 
MSTO will decrease, so for a fixed RGB mass loss rate, 
the older (i.e., less massive) stars will have smaller post- 
RGB envelope masses and thus will be hotter HBs. As 
a result, the UVX has been proposed as a p ossible age 
indicator for the centers of elliptical galaxies (|Greggio fe| 
Renzini|1999||Bressanet al.|1994||Yi et al |1999| ). |Bressan| 

LiB-mar 



et al 

power the 
than ~ 7.6 x 



(|1994| Lund that the'EHB'and AGi 

UV flux of a high mass galaxy at ages older 



manque stars 



the M31 bulge could therefore potentially explain the 
gradient in HP-HB stars we measure. 

There is, however, no o bvious evidence for an age gra- 
dient in the bulge of M31. Sagliaetal. (2010) constrained 
the stellar populations in the inner bulge of M31 using 
spline-interpolat ed Lick indices mode ls with a-element 
overabundances ( Thomas et al. 20031), assuming simple 
stellar populat ion models (Maraston 1998, 2005) and a 
Kroupa] ( |2001[ ) IMF. They found the age oFThe bulge 
(over the regions we analyze) to be consistent with the 
age of the universe, with no sig nificant radial grad ient. 
If anything, the mean ages from |Saglia et al. ( 2010[ ) sug- 
gest a positive age gradient with radius (i.e., older stars 
at larger radii) , which is in the opposite sense of what is 
needed to produce the observed HP-HB gradient. 

Another age effect is expected if the UVX is not from 
EHB stars and their progeny at all, but is instead from 
binary stars. This is an attractive explanation for why 
the surface density profile of LMXBs matches that of 
the HP-HB stars in Figure | 9| Studies by Ha n and oth- 

^20071 and refer- 

TTTr 



ers (Han et al. 2002 2003 Han et al. 



ences therein) present detailed synthetic SEDs that in- 
clude three binary star evolutionary channels: common 
envelope, Roche lobe overflow, and merging He white 
dwarfs. Over time, stars in these models will power the 
UVX. However, if this is the case, all old p opulations 



shoul d have a UVX, which they do not (e.g., O'Connell 
1999|. Alternatively, the lack of a central peak ol PJNe 



may suggest an alternate channel of P-AGB evolution in 
the central regions, such as AGB-manquc or PE-AGB, 
which in turn produces a central peak in the HP-HB 
stars. These conclusions should be considered tentative 



in ligh t of the likely incompleteness of the Merrett et al. 
(12006} PNe catalog in this region of M31. 



4.2. Galactic Metallicity Gradient 

Negative radial metallicity gradients are co mmon fea- 
tures of ellipticals and large spiral bulges (e.g.,|Ro ediger 
et al.||20111 Uablonka et al ; lj2007l IPapovich et al.l|2001t 
Kobayashi fc Arimoto|19 99 Dav ies et al.|1993| with typ- 
ical values ol A[Fe/H|/Alog(r) ~ —0.5 lor ellipticals. 
Metal-rich stars are linked to more mass loss on the RGB 
( Greggio fc Renzini|1990 ) which result in smaller HB en- 
velopes and thus hotter stars and more UV flux. Thus 
metallicity gradients may influence the fraction of stars 
that enter the EHB phase at each radius. 

Support for the importance of RGB ma ss loss in pro- 
ducing the UVX population comes from |Kalirai et al.| 



10 9 yr. A negative age gradient across 



(2007), who found that the white dwarfs of the metal- 
richold open cluster NGC 6791 ([Fe/H]~ +0.4, t = 8 
Gyr) are under-massive due to the enhanced mass loss 
of their progenitors. The same cluster hosts a significant 
population of very hot helium burning stars, which, al- 
though only 30% of al l HB stars, would be sufficient to 
produce a UV exces s. Dorman et al.| ( 1995[ ); Buzzoni & 
Gonzalez-Lopezlira (2008) found with synthetic popula- 
tion models that only ~ 20% of HB stars need to be hot 
to explain the UV upturn, in agreement with the NGC 
6791 observations. 

The small fraction of evolving MS stars that go though 
the HP-HB channel could be explained by a galactic 
metallicity gradient if the HP-HB are indeed the result of 
ext reme mass loss, and the mass loss depends on metallic- 
ity. Saglia et al. (20101 found that the metallicity in the 
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Fig. 11. — Strong anti-correlation as a function of radius between 



GALEX color (blue; |Thilker et al,|2005|l and the Lick Index -based 
metallicity and dispersion derived by paglia et al.| ( |2010[ ) (right 



axis, black). There is a sharp increase in metallicity where the 
FUV — NUV color becomes extremely blue. This anti-correlation 
may be evidence of the galactic metallicity gradient driving the 
radial variation in the numbers of UV-bright stars detected. The 
extreme blue color further shows that the UV-bright populations 
are strongly confined to the inner ~ 300 pc (~ 75") of the M31 
bulge. 

bulge of M31 decreases by 0.2 dex every decade in radius, 
from [Z/Hj ~ 0.4 in the center to solar at ~ 100" (see 
Figure 10). In the inner regions of M31, it is reasonable 



to assume that the bulk of the stars will have metal- 
licities that sample the metallicities reported by Saglia 
et al.| (2010), but that there will also be a tail to high 
metallicities. If metallicity variations drive the gradient 
in HP-HB stars, then the HP-HB population must result 
from this high metallicity tail, to avoid more than 
of stars passing through the HP-HB channel (i.e., Figure 

To}. 

Ancillary support for metallicity gradient being the 
driver of UV-bright star production can be found in Fig- 
ure [ill where we plot GALEX FUV - N UV vs. ra- 
diusTa long with the metallicities derived by |Saglia et al.| 
d2010h. The FUV — NUV color should track the relative 



st rength of the EHB (se e GALEX globular cluster CMDs 
in Schiavon et al.|[2~012 ). It is clear that there is a sharp 
upturn in [Z/HJ exactly where the FUV — NUV colors 
become extremely blue. This anti-correlation suggests 
that the radial variation in the numbers of UV-bright 
stars follows the galactic metallicity gradient. 

An alternative metallicity-dependent explanation for 
the UV excess is through production of low metallicity 
(Z < 0.008) MSTO stars of intermediate age, which emit 
a large amount of their flux in the NUV. If these stars 
were responsible for the UV flux we detect, we would 
expect the UV flux to fall toward the center, due to the 
increased metallicity. In fact, we see the opposite gradi- 
ent, the UV flux increases toward the center (see B98 
for a similar discussion). This mechanism does not help 
to explain the changing fraction of HP-HB stars. 



4.3. Elemental Abundances Gradients 
Carter et al. (20111) found that the UVX depends more 



in [a /Fe] could be expected if high amounts of star for- 
mation and supernovae increase the q -element enrich- 
ment of s tars in the bulges of galaxies (i Gre ggio fe Ren-| 
zmi||1999[ |Bressan et al.||1994| |Catelan||2009p . However, 
Saglia et al.J ( |20lJ ) found no obvious gradient in [a/Fe] 
over the regions we analyze in M31. Within the bulge of 
M31, the negative radial metallicity gradient exhibits a 
much stronger correlation to the stellar populations re- 
sponsible for the UVX. 

Higher metallici ty also implies higher h elium c ontent 
(AY/AZ ~ 2-3 , |Greggio fe Renzinij 19901 2.5-5 , |Bres~ 
|san et al.||1994[ at least 2.5, |Pagel et al.||l992|). TEe 



|san et al.]|1994| at least 2.5, |Pagel et al.| 
MSTO mass decreases strongl y with increasin g helium 



(sec 



0'ConneU||T999l and ref- 



abundance for a given age 

erences therein). For a fixed amount of KGB mass loss 
increased helium content will lead to hotter HB stars, 
that is, more EHB stars. Higher initial abundances of 
helium will also cause the star to burn more of its hy- 
drogen envelope during the HB phase, also favoring the 



produ ction of HP-HB stars ( |0'Connellil999||Horch et al 
1992|. The increased helium abundances necessary to 



produce EHB stars could be from contami nation by th e 



winds of older generations of AGB stars (Norri s 2004) 
the winds of massive stars ( Decressin et al. 20U7p . Alter 



strongly on [a/FeJ than [Fe/H]. An abundance gradient 



natively, these stars may have formed f rom the ac cretion 
of material from winds of central stars ( Seth|20T0 ). High 
helium content is a necessary ingredient to reproduce our 
results, and stellar winds causing the higher helium could 
explain the radial gradients we detect. Detailed stellar 
evolution modeling that varies the possible range of he- 
lium abundances from wind contamination is necessary 
to conclusively explain the radial gradient in MS stars 
that become HP-HB stars. 

4.4. Radial Dust Gradients 

Differential dust reddening could lead to gradients in 
the number of UV-bright stars detected as a function of 
radius. It is evident from Figure [2] that there is some 
structured dust obscuration in the oulge. The question 
then becomes, is there systematically less dust in the 
center of the bulge and if so, could dust alone account 
for the observed increase in faint UV sources towards the 
center of the galaxy? 

Based on three lines of evidence, we find it unlikely 
that the UV population gradients are caused by dust. 
First, M31's bulge is a 3-dimensional object, and its dust 
is likely to be confined primarily to the midplanc. It is 
therefore reasonable to assume ~ 50% of the bugle stars 
are in front of the dust layer. Second, we have inspected 
far-IR images of the bulge and see no systematic gradient 
in the FIR surface brightness over the central kiloparsec. 
We also see no obvious dust h ole in the central regions. 
Finally, an A b extinction map ( Melchior et al.|2000 Cia- 
rdullo et al.||1988 ) shows no evidence of a significant ra- 
dial dust gradient in the inner regions of M31. Therefore, 
while it is likely that dust has reduced the total number 
of UV stars we have detected (such that the percentages 
in Figure 10 are likely to be lower limits), it is unlikely 
that radial variations in extinction can be responsible for 
the rapid decline in HP-HB stars with radius. 

5. CLOSING REMARKS 

We have shown that the UV-bright population in the 
bulge of M31 is due largely to the progeny of hot HB 
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stars. Our findings are consistent with the conclusion portant in making AGB-manque and PE-AGB stars. 



of |B98l that these stars are just the tip of the iceberg 
of the hot HB stars responsible for the UVX. We report 
radial gradients in the resolved UV stellar population of 
M31's bulge that may be imperceptible in other data 
(e.g., in optical/NIR observations) because they regard 
less than ~ 3% of the evolving stellar populations in 
the bulge (~ 2% for M31's inner arcseconds and < 5% 



for M32 cf., |B98| |Brown et al. 2000J ). We have demon- 
strated the great utility of PHAT wide-area data, though 
bluer, deeper data are necessary for directly measuring 
the UVX. 

The gradients in the likely HP-HB stars can only re- 
sult from gradients in the properties of their progenitors 
across the bulge. The most likely changes in progenitor 
properties across the bulge are due to gradients in metal- 
licity, possibly coupled to gradients in helium abundance. 
Detailed analysis of the most probable combinations of 
these parameters are in a forthcoming paper (Rosenfield 
et al., in prep.) where we use stellar models and pop- 
ulation synthesis to reproduce the ratio of P-AGB to 
AGB-manque stars. This will give more precise estimates 
on which stellar evolution parameters are the most im- 



This work was supported by the Space Telescope Sci- 
ence Institute through GO-12055. PR wishes to express 
his gratitude to L.G., A.B., and Paola Marigo for invalu- 
able guidance and for graciously hosting him in Padova. 
We thank Brent Groves and the M31 Herschel imaging 
collaboration for allowing us to examine early Herschel 
observations of the bulge. We are grateful to Pauline 
Barmby for providing us the Spitzer IRAC images. We 
thank Roberto Saglia for providing us with metallicity 
values of the M31 bulge. This research has made use of 
the NASA/IPAC Extragalactic Database (NED), which 
is operated by JPL/Caltech, under contract with NASA. 
The FOC data have been obtained from the Multimis- 
sion Archive at the Space Telescope Science Institute 
(MAST). STScI is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under NASA 
contract NAS5-26555. L.G. and A.B. acknowledge sup- 
port from contract ASI-INAF 1/009/10/0. 

Facilities: HST (UVIS). 



REFERENCES 



Barmby, P., et al. 2006, ApJ, 650, L45 

Bell, E. F., & de Jong, R. S. 2001, ApJ, 550, 212 

Bensby, T., et al. 2010, A&A, 512, A41+ 

Bertola, F., Bressan, A., Burstein, D., Buson, L. M., Chiosi, O, 

& di Serego Alighieri, S. 1995, ApJ, 438, 680 
Bressan, A., Chiosi, C, & Fagotto, F. 1994, ApJS, 94, 63 
Bressan, A., Marigo, P., Girardi, L., Salasnich, B., Rubele, S., & 

Nanni, A. 2012, A&A, submitted 
Brown, T. M., Bowers, C. W., Kimble, R. A., Sweigart, A. V., & 

Ferguson, H. C. 2000, ApJ, 532, 308 
Brown, T. M., Ferguson, H. C, Davidsen, A. F., & Dorman, B. 

1997, ApJ, 482, 685 
Brown, T. M., Ferguson, H. C, Stanford, S. A., & Deharveng, 

J.-M., B. 1998, ApJ, 504, 113 
Brown, T. M., Smith, E., Ferguson, H. C, Sweigart, A. V., 

Kimble, R. a., & Bowers, C. W. 2008, ApJ, 682, 319 
Burstein, D., Bertola, F., Buson, L. M., Faber, S. M., & Lauer, 

T. R. 1988, ApJ, 328, 440 
Buzzoni, A., & Gonzalez-Lopezlira, R. A. 2008, ApJ, 686, 1007 
Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 

245 

Carollo, C. M., Danziger, I. J., & Buson, L. 1993, MNRAS 
Carter, D., Pass, S., Kennedy, J., Karick, A. M., & Smith, R. J. 

2011, MNRAS, 656 
Catelan, M. 2009, Horizontal Branch Stars and the Ultraviolet 

Universe, Astrophysics and Space Science Proceedings 

(Springer New York), 1-16 
Ciardullo, R., Rubin, V. C, Ford, Jr., W. K., Jacoby, G. H., & 

Ford, H. C. 1988, AJ, 95, 438 
Code, A. D. 1969, PASP, 81, 475 

Courteau, S., Widrow, L. M., McDonald, M., Guhathakurta, P., 
Gilbert, K. M., Zhu, Y., Beaton, R. L., & Majewski, S. R. 
2011, ArXiv e-prints 

Dalcanton, J. J., et al. 2012, ApJS, 200, 18 

Davies, R. L., Sadler, E. M., & Peletier, R. F. 1993, MNRAS, 
262, 650 

Decressin, T., Meynet, G., Charbonnel, C, Prantzos, N., & 

Ekstrom, S. 2007, A&A, 1044, 1029 
Dolphin, A. E. 2000, PASP, 112, 1383 

Dorman, B., O'Connell, R. W., & Rood, R. T. 1995, ApJ, 442, 
105 

Ferguson, H. C, & Davidsen, A. F. 1993, ApJ, 408, 92 
Girardi, L., Groenewegen, M. A. T., Hatziminaoglou, E., & da 

Costa, L. 2005, A&A, 436, 895 
Girardi, L., et al. 2008, PASP, 120, 583 
Girardi, L., et al. 2010, ApJ, 724, 1030 



Gorgas, J., Efstathiou, G., & Aragon Salamanca, A. 1990, 

MNRAS, 245, 217 
Greggio, L., & Renzini, A. 1990, ApJ, 364, 35 
— . 1999, Mem. Soc. Astron. Italiana, 70, 691 
Han, Z., Podsiadlowski, P., & Lynas-Gray, A. E. 2007, MNRAS, 

380, 1098 

Han, Z., Podsiadlowski, P., Maxted, P. F. L., & Marsh, T. R. 

2003, MNRAS, 341, 669 
Han, Z., Podsiadlowski, P., Maxted, P. F. L., Marsh, T. R., & 

Ivanova, N. 2002, MNRAS, 336, 449 
Horch, E., Demarque, P., & Pinsonneault, M. 1992, ApJ, 388, L53 
Jablonka, P., Gorgas, J., & Goudfrooij, P. 2007, A&A, 474, 763 
Kalirai, J. S., Bergeron, P., Hansen, B. M. S., Kelson, D. D., 

Reitzel, D. B., Rich, R. M., & Richer, H. B. 2007, ApJ, 671, 748 
Kang, Y., Bianchi, L., & Rey, S.-C. 2009, ApJ, 703, 614 
King, I. R., Stanford, S. A., & Crane, P. 1995, AJ, 109, 164 
King, I. R., et al. 1992, ApJ, 397, L35 
Kobayashi, C, & Arimoto, N. 1999, ApJ, 527, 573 
Koekemoer, A. M., Fruchter, A. S., Hook, R. N., & Hack, W. 
2002, in The 2002 HST Calibration Workshop : Hubble after 
the Installation of the ACS and the NICMOS Cooling System, 
ed. S. Arribas, A. Koekemoer, & B. Whitmore, 337 — h 
Kroupa, P. 2001, MNRAS, 322, 231 
Maraston, C. 1998, MNRAS, 300, 872 
— . 2005, MNRAS, 362, 799 

McConnachie, A. W., Irwin, M. J., Ferguson, A. M. N., Ibata, 
R. A., Lewis, G. F., & Tanvir, N. 2005, MNRAS, 356, 979 

Melchior, A.-L., Viallefond, F., Guelin, M., & Neininger, N. 2000, 
MNRAS, 312, L29 

Merrett, H. R., et al. 2006, MNRAS, 142, 120 

Norris, J. E. 2004, ApJ, 612, L25 

O'Connell, R. W. 1999, ARA&A, 37, 603 

O'Donnell, J. E. 1994, ApJ, 422, 158 

Oh, S.-H., de Blok, W. J. G., Walter, F., Brinks, E., & Kennicutt, 

Jr., R. C. 2008, AJ, 136, 2761 
Pagel, B. E. J., Simonson, E. A., Terlevich, R. J., & Edmunds, 

M. G. 1992, MNRAS, 255, 325 
Papovich, C, Dickinson, M., & Ferguson, H. C. 2001, ApJ, 559, 

620 

Peletier, R. F., Davies, R. L., Illingworth, G. D., Davis, L. E., & 

Cawson, M. 1990, The Astronomical Journal, 100, 1091 
Renzini, A., & Buzzoni, A. 1986, in Astrophysics and Space 

Science Library, Vol. 122, Spectral Evolution of Galaxies, ed. 

C. Chiosi & A. Renzini, 195-231 
Rocdiger, J. C, Courteau, S., MacArthur, L. A., & McDonald, 

M. 2011, MNRAS, 1261 



Bright UV Stars in the Bulge of M31 



13 



Saglia, R. P., et al. 2010, A&A, 509, A61 
Schiavon, R. P., et al. 2012, AJ, 143, 121 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 
525 

Seth, A. C. 2010, ApJ, 725, 670 
Skrutskie, M. F., et al. 2006, AJ, 131, 1163 
Thilkcr, D. A., ct al. 2005, ApJ, 619, L67 

Thomas, D., Maraston, C, & Bender, R. 2003, MNRAS, 339, 897 
Thomsen, B., & Baum, W. A. 1987, ApJ, 315, 460 
Trager, S. C, Faber, S. M., Worthey, G., & Gonzalez, J. J. 2000, 
The Astronomical Journal, 119, 1645 



Vassiliadis, E., & Wood, P. R. 1994, ApJS, 92, 125 
Voss, R., & Gilfanov, M. 2007a, A&A, 468, 49 
— . 2007b, MNRAS, 380, 1685 

Walterbos, R. A. M., & Kennicutt, Jr., R. C. 1987, A&AS, 69, 311 
Wirth, A. 1981, AJ, 86, 981 

Yi, S., Lee, Y.-w., Woo, J.-h., Park, J.-h., & Oemler, A. 1999, 
ApJ, 1, 128 



